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Abstract 

The morphogenetic transition of motile coral larvae into sessile primary polyps is triggered and genetically programmed 
upon exposure to environmental biomaterials, such as crustose coralline algae (CCA) and bacterial biofilms. Although the 
specific chemical cues that trigger coral larval morphogenesis are poorly understood there is much more information 
available on the genes that play a role in this early life phase. Putative chemical cues from natural biomaterials yielded 
defined chemical samples that triggered different morphogenetic outcomes: an extract derived from a CCA-associated 
Pseudoalteromonas bacterium that induced metamorphosis, characterized by non-attached metamorphosed juveniles; and 
two fractions of the CCA Hydrolithon onkodes (Heydrich) that induced settlement, characterized by attached 
metamorphosed juveniles. In an effort to distinguish the genes involved in these two morphogenetic transitions, 
competent larvae of the coral Acropora millepora were exposed to these predictable cues and the expression profiles of 47 
coral genes of interest (GOI) were investigated after only 1 hour of exposure using multiplex RT-qPCR. Thirty-two GOI were 
differentially expressed, indicating a putative role during the early regulation of morphogenesis. The most striking 
differences were observed for immunity-related genes, hypothesized to be involved in cell recognition and adhesion, and 
for fluorescent protein genes. Principal component analysis of gene expression profiles resulted in separation between the 
different morphogenetic cues and exposure times, and not only identified those genes involved in the early response but 
also those which influenced downstream biological changes leading to larval metamorphosis or settlement. 
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Introduction 

Among the cues that trigger the pelago-benthic transition of 
coral larvae, biomolecular cues associated with the natural habitat, 
such as coral rubble, crustose coralline algae (CCA) and marine 
bacteria, have received considerable attention in the literature [1- 
4], Previous investigations of gene expression profiles in coral 
larvae after exposure to these natural morphogenetic cues have 
focussed on the regulation of candidate genes with presumptive 
key functions in coral development, such as those implicated in cell 
proliferation, apoptosis, differentiation, migration, adhesion, 
biomineralization and immunity [5—11]. These investigations 
identified a large number of genes that were differentially 
expressed across the complete developmental transition, from 
competent swimming larvae to post-metamorphosis primary 
polyps [7], and that regulatory changes in gene expression 
occurred between 4 to 12 h post incubation (hpi) to a variety of 
morphogenetic stimuli [6-8, 1 2] . The general conclusion was that, 
at the transcriptional level, coral larvae appear to anticipate 
metamorphosis [6]. 



Changes in gene expression profiles of larvae exposed to 
morphogenetic stimuli have mostly been investigated using non- 
targeted approaches [6,8] primarily due to the variable and 
complex composition of coral rubble obtained in the field, 
typically characterized by a mix of live and dead CCA, endolithic 
green algae, polychaetes, diatoms, protozoa and bacteria. Further 
compounding this issue is the likelihood that changes to the gene 
expression profiles of larvae exposed to CCA are in fact an 
integrated response to a multitude of qualitatively different triggers 
possibly affecting gene expression at multiple and unrelated levels. 

A bioassay-guided chemical isolation strategy was employed to 
generate chemically refined and quantifiable putative morphoge- 
netic cues from these natural biomaterials that would trigger larval 
morphogenesis in vitro. This strategy yielded (1) an organic and an 
aqueous fraction of the CCA Hydrolithon onkodes (Heydrich) that, 
under laboratory conditions, triggered in a highly reproducible 
manner, settlement, characterized by attached metamorphosed 
juveniles [13] and (2) a crude extract of Pseudoalteromonas strain 
JO 10 (J010-E) isolated from H. onkodes, known to contain 
tetrabromopyrrole (TBP), that triggered larval metamorphosis, 
characterized by non-attached metamorphosed juveniles [4, 1 2] . 
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The conceptual approach taken in the present study targeted 
genes of interest (GOI) putatively involved in the early stages of 
coral larval morphogenesis (1-3 hpi), using multiplex RT-qPCR 
[12] to measure the gene response after short-term exposure to the 
chemically refined morphogenetic cues, and ultimately to 
discriminate between their possible modes of action. This 
approach had several advantages over previous studies: (1) it 
reduced the variance within biological replicates [14], (2) 
improved the resolution of gene expression levels, (3) lowered 
possible interferences from other components normally present in 
environmental biomaterials (see above) that may elicit non-specific 
changes in gene expression levels [7], and (4) decreased secondary 
time-related effects of cues on downstream gene expression. 
Discussed here are the putative roles of GOI in the early stages of 
coral larval settlement and metamorphosis based on changes in 
their expression profiles prior to any observed morphogenetic 
transformation. 

Materials and Methods 

Ethics Statement 

This work was undertaken under permit G 10/33440.1, issued 
by the Great Barrier Reef Marine Park Authority. No ethical 
approval was required for any of the experimental research 
described herein. 

Sample collections and larval maintenance 

Colonies of the scleractinian coral A. millepora were collected 
from Pelorus Island on the Great Barrier Reef (GBR; 18°33'S 
146°29'E) and transported to the Australian Institute of Marine 
Science (AIMS) prior to the predicted spawning events in 
November 2010 and November 2011. Chips of the CCA 
Hydrolithon onkodes (Heydrich) were collected from Trunk Reef 
(GBR; 18°24'S, 146°48'E) and transported to AIMS in November 
of 2010, 2011 and 2012. The taxonomic identity of aU CCA 
samples was confirmed prior to use [15]. Each CCA chip was 
carefully removed from the solid substrate using sterile bone 
cutters creating a thin layer without any other visible algae or 
invertebrates and subsequently washed with 0.2 um filtered 
seawater (FSW). Coral larvae were obtained after fertilization of 
gametes as described previously [4]. Larvae were raised in flow- 
through 500 L tanks containing FSW at 27-28°C until they 
became fully competent, as characterised by at least 80% of larvae 
displaying an elongated shape, exhibiting active benthic searching 
behaviour and settlement upon exposure to CCA. In this 
controlled environment larvae are deemed competent approxi- 
mately four days after fertilisation and can remain in this state for 
an additional ten days [12]. Their high lipid content is used for 
buoyancy and as an energy source [16] allowing them to survive at 
least three weeks after settlement without feeding. 

Growth conditions and extraction of metamorphosis- 
inducing bacteria 

The larval metamorphosis-inducing cue was isolated from 
Pseudolalteromonas strain JO 10 (JF309049) [4]. Briefly, bacterial 
colonies (3.2 g) were carefully scraped off agar plates, suspended in 
50 mL ethanol (64 mg/mL) and sonicated for 20 min (Soniclean 
500T, SA). The resulting bacterial extract (J010-E) was filtered 
(0.2 um) and stored in the dark at -20°C. 

Preparation of CCA-derived settlement cues 

Two settlement-inducing cues, derived from an organic extract 
and an aqueous extract, were isolated from H. onkodes. An ethanol 
extract of H. onkodes (2300 g wet weight, November 2011) was 



fractionated on reversed phase C-18 vacuum flash chromatogra- 
phy (RP-18 VFC, 80x600 mm) by elution with Milli-Q water 
(MOW, 600 mL), 10% methanol in MQW (600 mL), 20% 
methanol in MQW (600 mL), 80% methanol (600 mL) in 
MOW and 100% methanol (1200 mL). The 80% methanol 
fraction (0.95 mg/mL; herein referred to as HO org ) induced high 
rates of larval settlement in preliminary bioassays and was stored 
at -20°C [13]. 

The aqueous extract was prepared as follows: H. onkodes (500 g 
wet weight, November 2011) was washed with 500 mL MQW 
three times to remove salts. These washings were discarded. The 
washed CCA were autoclaved at 121°C for 60 min in 250 mL 
MQW. This aqueous extract was decanted and the algae 
autoclaved twice more (121°C for 60 min in 250 mL). The 
pooled extract was ultrafiltrated (Grace Microcon-100, > 100 kD), 
freeze-dried, resuspended in 250 mL 10% ethanol (56.1 mg/mL) 
in MQW, herein referred to as HO aq , and stored in the dark at 
4°C [13]. 

HO ()rg was extracted over three consecutive years (2010-2012); 
HO aq was extracted over two consecutive years (2011-2012). 
Comparison of the proton nuclear magnetic resonance spectra of 
each of these samples confirmed their chemical composition was 
consistent (data not shown). 

Determination of minimum exposure time of larvae to 
cues that elicit settlement and/or metamorphosis 

Preliminary experiments were conducted using 6-9 day-old 
swimming larvae in 6-well plates. In 2010, larvae were exposed to 
J010-E (Table 1) to assess the minimum required exposure time of 
larvae to the metamorphosis-inducing cue. The extract J010-E 
(2 uL; 128 ug/mL) or 100% ethanol (2 uL negative control) were 
added to 10 mL of 0.2 um FSW, mixed thoroughly and 
transferred to 6-well plates (Nunclon Delta Surface, Thermo 
Scientific, Denmark; final concentration 12.8 Ug/mL). Larvae 
were added to each well (n = 20 per well) and monitored hourly for 
visible signs of metamorphosis without settlement, characterised by 
flattening into discs and septal development in floating polyps. 
Following the same protocol, swimming larvae were exposed to 
settlement-inducing cues HO org (evaporated on the dish surface 
before adding FSW) or HO aq (aliquot added to FSW) in 201 1. No 
cue was added to the control treatments (n = 6). In a preliminary 
experiment, addition to the well and evaporation of the carrier 
solvents (10% ethanol, 100% ethanol, 80% methanol and MQW) 
to dryness resulted in no settlement, as observed with the no 
solvent (blank) control (A. Negri, personal communication), and 
therefore only the blank control was used. 

Experimental design and sample collection 

Gene expression profiles may vary considerably throughout the 
different coral life phases; pre-settlement motile larvae to post- 
settlement metamorphosed primary polyps [7]. To avoid any 
variation due to differences in larval development and to 
accomplish meaningful comparisons between larvae exposed to 
control conditions versus J010-E or CCA-derived HO ()rg and 
HO aq cues, all assays were terminated before the onset of visible 
metamorphosis without attachment or settlement (as determined 
above). Assays with J010-E were terminated 1, 2 and 3 hours post 
incubation (hpi), whereas assays with HO org and HO aq were 
terminated after 1 hpi (see experimental parameters in Table 1). 

The extract J010-E (60 uL; 3.84 mg/mL) or 100% ethanol 
(60 uL, control) was mixed with 250 mL FSW in 400 mL glass 
beakers (Boeco, Germany; 6 replicates per treatment) to ensure 
homogeneous distribution prior to the addition of a further 50 mL 
FSW containing competent swimming larvae (six days-old, 
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Table 1. Experimental parameters. 





Cue 


HO 

""org 


HO 


J010-E 


Origin of Cue 


Hydrolithon onkodes (Heydric 


h) crustose coralline alga 


Biofilm of Pseudoalteromonas strain J010 




(CCA) 




(JF309049) 


Larval response 


Settlement & metamorphosis 




Metamorphosis without attachment 


Spawning Year 


2011 




2010 


Larval age 


6 days 




9 days 


Larval numbers 


100-200 




200-300 


Treatment volume 


100 mL 




300 mL 


Exposure time 


1 hpi 




1-3 hpi 


Replicates 


6 




18 (6 per h) 


Aliquot (mD 


25 |iL 




60 uL 


Application 


Dry 


FSW 


ethanol 


Final cue concentration 


238 ng/mL 


14 (ig/mL 


12.8 (ig/mL 


Settlement/metamorphosis rate 


16 hpi 94±5% 


16 hpi 70±8% 


6 hpi 98 ±2% 


Assay 1 - Most stable genes 


Rps7, ATF4/5 




Rpl9, CTL-2 


Average %CV ± STDEV (assay 1 ) 


6.9±5.1 




8.2±5.5 


Assay 2 - Most stable genes 


Rps7, Amgalaxin-like-1 




amilRFP, Rps7 


Average %CV ± STDEV (assay 2) 


9.6±5.8 




7.7±6.5 



doi:1 0.1 371 /journal.pone.0091 082.t001 



n~300, final concentration 12.8 ug/mL). Aliquots of HO org 
(25 uL; 23.8 ug/mL) were applied to the base of 400 mL glass 
beakers and allowed to evaporate, while HO aq (25 uL; 1.4 mg/ 
mL) was applied directly into the seawater (6 replicates per 
treatment). FSW (50 mL) was added to 400 mL glass beakers and 
gently mixed prior to the addition of 50 mL of FSW containing 
swimming larvae (nine days-old, n~200) to give a final volume of 
100 mL (final concentration HO org 238 ng/mL and HO aq 14 [ig/ 
mL). Blank treatments (no cues or solvents) were used as the 
negative controls. 

In all experiments, an additional glass beaker containing each of 
the cues was left for 24 hpi to monitor larval metamorphosis and 
mortality. After each observation time point, larvae were collected 
on 5 |im, 25 mm diameter sterile filters under vacuum. Retained 
larvae were transferred into 1.5 mL cryovials, snap-frozen in 
liquid nitrogen and stored at — 80°C. 

Multiplex RT-qPCR Assays 

Forty-seven genes of interest (GOI) were selected for their 
putative role(s) in larval settiement or in coral immune defence 
pathways, and incorporated in two multiplex RT-qPCR assays 
(Assay 1; immunity-related genes potentially involved in ontogen- 
esis, and Assay 2; genes reportedly involved in settlement; 
kanamycin (Kan) internal control gene in both assays) as previously 
described (Tables SI and S2 in [12]). The suggested roles of the 
GOI in relation to marine invertebrate larval settiement and 
metamorphosis have been summarised and reported previously 
(Table 1 in [12]). One negative template (T— ) and one negative 
Reverse Transcriptase (RT— ; Kan only) control was used in each 
assay to test for cDNA and RNA contaminants, respectively, 
according to the GenomeLab™ GeXP manual (a detailed 
description for this procedure is reported in [17]). For technical 
reasons the number of genes in every multiplex assay is limited to 
approximately 30 [12,17]. The protocols for extraction of mRNA, 
preparation of cDNA, amplification by PCR and subsequent 
electrophoresis were adopted from Siboni et al. [12] with 
modifications: 1) larval mRNA was measured twice (NanoDrop 



ND-1000, Thermo Scientific, USA), averaged and diluted to 
5 ng/uL (according to the GenomeLab™ GeXP manual); 2) 
10 ng of mRNA was used in the reverse transcription mix; and 3) 
dilution of PCR products was changed to 1:40 forJOlO-E and 1:20 
for HO cues. Diluted PCR products were analysed individually on 
an automated capillary electrophoresis sequencer (CEQ™ 8800 
Genetic Analysis System, Beckman-Coulter, Fullerton, CA, USA) 
as follows: diluted PCR products (1 mL) were loaded into a 96- 
well plate with 0.2 uL of DNA size standard (400 bp; Beckman- 
Coulter) and 38.8 uL of Genome-Lab Sample loading solution. 
Electropherograms were visualized, filtered and matrices generat- 
ed for further analyses following the GenomeLab™ GeXP 
manual (a detailed description for this procedure is reported in 
Souter et al. [17]). 

Statistical Analysis 

Data were normalised against Kan . An expression stability 
measure for all genes in both assays was performed using the 
GeNorm program (http://medgen.ugent.be/genorm/) according 
to Vandesompele et al. [18]. The most stable pair of genes was 
used to normalize the gene expression levels of the GOI in 
Microsoft Excel (2007) using the Geomean function (Table 1). 
Further interrogation of the data was undertaken by normalisation 
using three of the most stable genes in every assay (Figure SI) and 
results compared with those from the pair of stable genes. For data 
analysis, only samples with at least two technical repeats and 
coefficient of variation percentage (%CV) lower than 25% were 
considered. Statistical analyses were conducted using Kruskal- 
Wallis (non-parametric one way analysis of variance by ranks) 
followed by multiple comparisons of mean ranks for all groups 
with STATISTICA version 10.0 (StatSoft, Tulsa, OK). Percent- 
age differences of those genes that were differentially regulated 
from control were calculated. Further, data were analysed by 
principal component analysis (PCA) with PAST [19]. All PCA's 
were performed on a variance-covariance matrix, as all variables 
were measured in the same units, and without any group 
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Change in gene regulation as % of control 



Gene 


Extracts of 
Hydrolithon Onkodes 


Bacterial TBP-containing extract (J010-E) 


Function 


Assay 1 


HOorglh 


HO aq 1 h 


1 h 


2h 


3h 


12 h 




TIR-1 










+42.9 


+43.2 


Toll/TLR pathway 


TRAF-6 










+22.1 


-35.6 


(MAPK and TIR pathways) 
alarm genes, cell cycle 


MEKK-I 




+79.8 


| +22.1 


+52.6 


+91.1 


ERK-2 




+32.9 






+ 15.6 




Regulation and proliferation 


MAPK-p38 




-73.3 




-21.8 


-22.4 




AP-l/cFos 


+58.1 


+71.8 


+43.6 


+44.7 


+62.9 


+ 14.8 


AP-1/c.Jun 


-64.6 


-45.2 


+45.8 


+82.4 


+ 129.3 


-21.7 




ATF4/5 








-11.1 






NF-kB 






-14.1 


-26.5 


-14.7 




TRAM | 


-37.4 










Apextrin 


-58.2 


-49.4 


+ 17.6 






-56.4 


MAC/PF domain, opsonization 
And cell recognition 


Tx60-A 


-58.3 


-41.5 






+41.7 




Bf 




-80. 1 










C3 


+46.2 


+29.5 




-15.7 






millectin 


-66.7 


-66.9 


+46.5 


+ 17.1 






Cell recognition and adhesion 


CTL-1 




-72.6 


+75.0 


+65.0 


+60.1 




CTL-2 


-17.5 










HL-2 


-79.6 


-83.7 


+ 182 


+66.0 


+53.5 


HL-3 


-79.9 


-74.1 


■7 9 


-13.6 






Assay 2 


Septin-J 








+ 12.4 






Cell division and cytokinesis 


Septin-7 


+19.0 




+29.5 


+27.4 




LRRC59 




+ 13.0 




+14.0 




+20.4 


Toll/TLR pathway 


CAMP/ATF2 


+27.9 






+ 12.4 


+ 14.4 




CAMP/ATF4 


+52 1 










-1 


cAMP'ATF6 










+66.7 


-36.4 


GABA/GAT-2 






+8.7 


+ 18.1 


+52.1 


+ 104 


Neurotransmitters 


GABA/GA T-3 






+ 18.2 




+41.8 


+ 120 


GABAunc-47 












+213 




GABA/ine 






+40.6 


+36.2 




NADH-ubiq 




+32.3 








Metabolism 


Lectin MemgalRl <■ 






+ 12.7 


+ 14.6 




Nematocyst development 


Lectin nemgalB ; 








+9.7 






Pha-C3 






+12.5 


+76.5 


+ 100 


+70.6 


cAMP-dependent protein kinase 


amilGFP \ i 


+ 11.1 




-16.7 


-23.0 


-35.8 


Fluorescent protein genes 


amilCFP 








-11.7 


-26. 1 


-64.2 


amilRFP 




+9.8 






i 3.4 


-36.5 


amilCP 




+29.0 








-21.4 


Colourless GFP-like 








-7.0 




-35.8 


Amgalaxin-like-I 




+6.8 






+ 176 


Taxonomically restricted genes 
calcium handling 


SCRiP-like 












-59.6 


SCRiP2 








-9.2 




-74.0 




SCRiP3 








+353 


+245 


+209 


Actin 




+33.8 








-40 2 


Motility 


Integrin-alphal | : \ 












Cell adhesion 



Figure 1. Differential gene expression following exposure to bacteria- and CCA-derived cues. Significant (p<0.05) change in gene 
regulation is given as the % difference compared to the control. Dark and light green (+) represents those genes that were up-regulated (X>20, 
20<X<0); orange and red (— ) represents those that were down-regulated (0>X>— 20, X<— 20). Data from the 12 hpi experiment (high 
concentration treatment and complete metamorphosis) was taken from Siboni et al. [12], which also includes full protein names and description of 
the genes. 

doi:10.1371/journal.pone.0091082.g001 
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assumptions. PCA score plots are presented with convex hulls 
highlighting groupings and showing PCI and PC2. 

Results 

Determination of minimum exposure times of larvae to 
cues to elicit settlement or metamorphosis without 
attachment 

The minimum exposure time to elicit larval responses varied 
among the different cues. Visible signs of early metamorphosis 
(flattening into discs and development forming floating polyps) 
occurred 4 hpi with J010-E, whereas complete metamorphosis 
without attachment of 98 ±2% of larvae occurred 6 hpi. In order 
to sample swimming larvae prior to any observable morphogenetic 
changes, they were collected 1-3 hpi with J010-E (Table 1). 

While no behavioural or developmental response was recorded 
during the first 2 hpi with HO org and HO aq , initial attachment was 
observed after 3 hpi, and complete metamorphosis with attach- 
ment was observed 16 hpi (settlement rates of 94±5% (HO org ) and 
70±8% (HO aq ), respectively). To reduce any variation in the 
assays and to ensure that all larvae were sampled at the same 
developmental phase, larvae were collected 1 hpi with HO org or 
HO aq (Table 1). It was noted that addition of the CCA-derived 
cues directly to larvae in FSW resulted in ~5% mortality, as 
determined by the release of mucus resulting from tissue 
breakdown. This method of application was deemed unsuitable 
for the current study. Pre-conditioning the FSW with the CCA- 
derived cues followed by addition of larvae did not lead to any 
mortality (as determined above by release of mucus); therefore 
these conditions were used in the main experiment. 

The larvae exposed to the three different cues displayed normal 
swimming behaviour with no morphological signs of setdement at 
the time of termination of experiments (1-3 hpi). With respect to 
the two CCA-derived cues, no larval response was recorded during 
the first 2 hpi. Furthermore, a subset of larvae setded, deposited a 
calcareous skeleton, divided and underwent complete metamor- 
phosis. These juveniles were successfully infected with Symbiodinium 
(data not shown). No mucus resulting from tissue breakdown was 
observed 24 hpi in response to any of the three cues. 

Gene regulation/expression 

Gene pairs used to normalize expression levels and their 
coefficients of variation (%CV) are listed in Table 1. For 11% of 
the total number of data points in this study two technical 
replicates were used instead of three, the third being unreliable. 
The negative template (T— ) control did not produce any 
measureable products, as described by the GenomeLab™ GeXP 
manual, while a clear measureable peak at 325 nucleotide size was 
observed for the negative reverse transcriptase control (RT— ; Kan 
only), confirming that electropherograms resulting from the larval 
samples reflected transcribed gene amplification products. Only 
genes that changed significandy [p< 0. 05, Kruskal-W allis followed 
by multiple comparisons of mean ranks) as compared to the 
control were further analysed (Figures 1 and SI). Two boxplot 
figures showing the variances, one for the CCA-derived cues, 
HO org and HO aq (Figure S2) and the other for the bacterial cue 
J010-E (Figure S3), represent the complete data. The averaged 
data is presented in Figure 1 . The inclusion of a third gene in the 
normalisation had only a minor impact on the results (Figure SI), 
therefore only the two most stable genes, were used (Figure 1). It 
should be noted that the original reference genes included in the 
two assays [12] were not always the most stable genes according to 
GeNorm program [18] (Table 1). In all instances at least one of 
these stable genes was a ribosomal protein gene. Given that 



significant differences in gene expression were observed at the 
early time point and that no mucus strands, indicative of tissue 
breakdown, were observed 24 hpi with any of the cues, the 
concentrations used in this study were considered appropriate to 
measure gene expression levels while maintaining larval compe- 
tency. 

Mutiplex RT-qPCR Assay 1; putative settlement/ 
metamorphosis-related immunity genes 

Expression levels of 16 GOI [MEKK-1, ERK-2, MAPK-p38, AP- 
1/cFos, AP-l/cJun, NF-kB, TRAM, Apextrin, Tx60-A, Bf, C3, 
milkctin, CTL-1, CTL-2, HL-2 and HL-3; Figure 1) changed 
significandy 1 hpi. During this time interval the lectins milkctin, 
CTL-1, CTL-2, HL-2 and HL-3 were predominately up-regulated 
in response to the metamorphosis-inducing cue JO 1 0-E, whereas 
their down-regulation was observed in response to one or both of 
the settiement-inducing CCA-derived cues HO org and HO aq . AP- 
1 IcFos was consistendy up-regulated in response to JO 1 0-E, HO ()rg 
and HO aq . The genes TIR-1, TRAF-6, MEKK-1, ERK-2, MAPK- 
p38, ATF4/5, Tx60-A and C3 (Figure 1) caused significant 
differences 2 and/or 3 hpi with JO 1 0-E, whereas AP I I cFos, AP- 
1 /cjun, JVF-kB, CTL-1 and HL-2 were significandy regulated at all 
three sampling time points (Figure 1). Some genes were similarly 
regulated following exposure to HO ()rg and HO aq (MAPK-p38, AP- 
1/cFos, AP-l/cJun, TRAM, Apextrin, Tx60-A, C3, milkctin, HL-2 and 
HL-3) while others were differentially regulated (as compared to 
control) after exposure to only one of the HO cues [MEKK-1, ERK- 

2, Bf, CTL-1 and CTL-2; Figure 1). The PCA score plot (Figure 2A) 
explained 79% of the variance and revealed clear differences in 
gene regulation following exposure to JO 1 0-E or the CCA-derived 
cues. Further PCA considering only the different exposure times to 
JO 1 0-E (Figure 2C) and another comparing the organic and 
aqueous CCA-derived cues (Figure 2E), explained 73% and 91% 
of the variance, respectively. There were clear differences between 
the time points, particularly 1-3 hpi compared to 12 hpi with 
JO 1 0-E, and between the two CCA-derived cues. The associated 
loading plots are available in (Figure S4). 

Mutiplex RT-qPCR Assay 2; other putative settlement/ 
metamorphosis-related genes 

Similar to Assay 1, expression levels of 16 GOI changed 
significandy after 1 hpi with at least one of the cues [Septin-7, 
LRRC59, cAMP/ATF2, cAMP/ATF4, GABA/GAT-2, GABA/GAT- 

3, MADH-ubiq, Lectin / nemgalRl , Pha-C3, amilGFP, amilRFP, amilCP, 
Amgalaxin-like- 1 , SCRIP3, Actin and Integrin-alphal; Figure 1). The 
fluorescent protein genes [amilGFP, amilRFP and amilCP) were 
significandy up-regulated in response to the CCA-derived cues 1 
hpi, and all FP genes were significantly down-regulated in response 
toJOlO-E but only after 2, 3 and/ or 12 hpi (Figure 1). Some genes 
presented significant differences only after 2 and/ or 3 hpi with 
J010-E (Septin-1, Septin-7, LRRC59, cAMP/ATF2, cAMP/ATF6, 
GABA/ine, Lectin /nemgalRl, Lectin I nemgalBl , amilGFP, amilCFP, 
amilRFP, Colorless GFP-like and SCRiP2; Figure 1), while others 
presented constant significant differences at all three sampling time 
points [GABA/GAT-2, Pha-C3, Amgalaxin-like- 1 and SCRIP3; 
Figure 1). Two genes [cAMP I ATF4 and amilGFP) were consistendy 
up-regulated following exposure to the two CCA-derived cues. 
Four genes [Septin-7, cAMP/ATF2, Lectin /nemgalRl and Integrin- 
alphal; Figure 1) were significantly up-regulated after exposure to 
only HO ()rg while six genes [LRRC59, NADH-ubiq, amilRFP, 
amilCP, Amgalaxin-like- 1 and Actin) were significantly up-regulated 
following exposure to HO aq . As for Assay 1, the PCA score plot 
(Figure 2B) explained 88% of the variance, and again revealed 
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Figure 2. Influence of different settlement-inducing cues and exposure times on gene expression of Acropora millepora larvae. PCA 
score plots of the % difference in gene expression compared to the control, presented with convex hulls highlighting groupings and showing PCI 
and PC2, for A) J010-E 1-12 hpi and H0 org/aq , Assay 1, B) J010-E 1-12 hpi and H0 org/aq , Assay 2, C) J010-E 1-12 hpi, Assay 1, D) J010-E 1-12 hpi, Assay 
2, E) HO org/aq , Assay 1 and F) HO org/aq , Assay 2. Time points for J010-E are represented by: 1 hpi = black full circle, 2 hpi = blue empty square, 3 
hpi = green empty circle and 12 hpi = cyan empty triangle. CCA-derived cues are represented by: HO org = pink full square and HO aq = red cross. 
doi:1 0.1 371 /journal.pone.0091 082.g002 



clear differences in gene regulation following exposure to JO 1 0-E 
or the CCA-derived cues. Additional PCA of the JO 1 0-E 
treatments only (Figure 2D) and of the CCA-derived cues only 
(Figure 2F) explained 84% and 70% of the variance, respectively. 
Again the PCA showed clear differences between the early time 
points (1-3 hpi) and 12 hpi [12], and some separation between the 
two CCA-derived cues. The associated loading plots are available 
in (Figure S5). 



Discussion 

Biofilms of Pseudoalteromonas strain JO 10, an extract of this 
bacterium JO 1 0-E) and the purified bacterial metamorphic 
metabolite, tetrabromopyrrole (TBP), all induce metamorphosis 
of A. millepora larvae without attachment [4]. In contrast, the CCA- 
derived cues HO org and HO aq , repeatedly extracted from H. 
onkodes over three consecutive years (2010-2012), reproducibly 
induced complete settlement. Substantial changes in expression 
profiles of 32 genes of interest (GOI; Figures 1 and 2) hypothesized 
to be involved in larval settlement were observed following one 
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hour post incubation (hpi) with each of these cues (J010-E, HO org 
or HO aq ) and prior to any visible signs of morphogenesis. The 
PCA score plots (Figure 2) show a clear separation between the 
different cues and time of exposure, thus correlating the early 
response of genes with downstream biological changes leading to 
metamorphosis or settlement. Statistical separation of the JO 1 0-E 
treatments from the GCA-derived cues HO ()rg and HO aq was 
largely explained by the regulation of the immunity-related genes 
MAPK-p38, API /cjun, the lectins millectin, HL-2, tissue remodel- 
ling genes Apextrin, Tx60-A, cysteine-rich protein genes SCRiP3, 
Amgalaxin-like-1 and Pha-C3 (Figures 2A-B, S4A and S5A). It 
should be noted that genes with a low expression difference as 
compared to the control, even though significant, will most likely 
have lower biological importance at that early time point in larval 
development. 

MAPK-p38 and AP-l/cJun form the AP-1 early response 
transcription factor of the Toll/TLR pathway [20] and are 
associated with the activating transcription factors cAAIP-ATF. All 
three cues showed a similar trend in the regulation of MEKK-1, 
ERK-2, MAPK-p38, AP-l/cFos, and cAMP/ATF2 (Figure 1), 
indicating common activation of the immunity response in the 
early developmental (1-3 hpi) phase [5,6,8,10-12,21,22]. Howev- 
er, AP-1 / cjun exhibited an opposite response to JO 1 0-E over time 
but showed the same response at 12 hpi and with the two CCA- 
derived cues. TRAM and NF-kB (Figure 1) were also differently 
expressed, with TRAM responding to both CCA-derived cues and 
NF-kB responding to JO 1 0-E, suggesting different roles during the 
early stages of settlement and metamorphosis. 

Changes in swimming and searching behaviour followed by 
larval settiement [1,6] involve cell-surface recognition mechanisms 
mediated by lectins [23] including substrate selection [8], 
activation of tissue remodelling [5,6,24], nematocyst development 
[5] and establishment of symbiosis [5]. Initial gene expression 
patterns in larvae exposed to the CCA-derived cues were 
characterised by down-regulation of the lectins millectin, CTL-1, 
CTL-2, HL-2 and HL-3 (Figure 1). Exposure to JO 1 0-E, however, 
resulted mostly in up-regulation (1-3 hpi) with gene expression 
returning to baseline levels after metamorphosis (12 hpi) 
highlighting opposing roles of these genes in the early stages of 
metamorphosis (1-3 hpi) and settiement. 

Strong up-regulation of Apextrin and Tx60-A during larval 
development correlates with tissue remodelling function 
[5,7,8,10,25]. Exposure to both CCA-derived cues 1 hpi resulted 
in down-regulation of both genes, indicating no initiation of tissue 
re-modelling. Exposure to JO 1 0-E elicited the opposite response 
suggesting tissue remodelling in preparation for metamorphosis 
had commenced. After 12 hpi with J010-E down-regulation of 
Apextrin was observed correlating with complete metamorphosis 
without attachment [12]. 

Galaxin-related proteins are important structural proteins in the 
calcifying organic matrix of scleractinian corals [26]. Amgalaxin- 
like-1 is expressed strongly in larvae and settled primary polyps but 
not in adult colonies [11]. The up-regulation of Amgalaxin-like-1 in 
larvae exposed to JO 1 0-E and HO aq (1 hpi to 12 hpi) confirms its 
involvement in early development. The expression of genes 
encoding for other cysteine-rich proteins, the SCRiPs proteins, 
also changed throughout the developmental transition suggesting 
distinct roles in coral development [22] and bio-mineralization 
[22]. SCRiP3 was up-regulated following 1-12 hpi with JO 1 0-E, 
while SCRiP2 was down-regulated in response to JO 1 0-E, similar to 
that reported previously [6-8,22] highlighting a role in metamor- 
phosis. Pha-C3, known to mediate larval development and 
settlement in Hydroides elegans [27] was up-regulated after exposure 
to JO 1 0-E. The expression levels of SCRiPs, Amgalaxin-like-1 and 



Pha-C3 genes did not return to baseline levels, even after 
metamorphosis in the water column, but continued to be regulated 
12 hpi with JO 1 0-E [12] correlating with continued development 
of septa. No change in Pha-Ci or SCRiPs expression profiles was 
observed after short term (1 hpi) exposure to either of the CCA- 
derived cues. 

Although Actin did not influence significantly the separation 
between JO 10-E and the CCA-derived cues, long-term (12 hpi) 
exposure to JO 1 0-E resulted in its down-regulation coinciding with 
the formation of non-swimming floating metamorphosed larvae 
[12], similar to that observed by Hayward et al. [7] while, up- 
regulation occurred in larvae exposed to HO aq (1 hpi). Further- 
more, the chemical composition of the two CCA-derived cues 
elicited different gene expression profiles 1 hpi even though both 
cues ultimately resulted in normal settlement 16 hpi. These results 
provide further evidence that up-regulation of Actin, which is 
linked to movement in coral larvae, likely drives larval settlement. 

Differential gene expression of swimming larvae exposed to 
JO 10-E over time (1-3 hpi and 12 hpi) enabled further decoupling 
of GOI involved in different stages of metamorphosis. In both 
assays, the PCA score plots associated with JO 10-E only 
(Figure 2C-D), showed that the early time points were clearly 
distinct from the 1 2 hpi time point and correlated with swimming 
and metamorphosed larva, respectively. This separation was 
strongly influenced by the regulation of several genes involved in 
the Toll/TLR pathway, cell recognition, cell adhesion, neuro- 
transmission, FP-related metabolism and calcium handling: AP-1 1 
cjun, MEKK-1, Apextrin, HL-2, GABA/GAT-3, Amgalaxin-like- 1 , 
SCRiP2, SCRiP3, with a lesser influence by T1R-1, TRAF-6, AP- 
1/cFos, millectin, CTL-2, Lectin/ nemgalRl , Lectin InemgalBl, amilGFP 
and amilCFP (Figures 2C-D, S4B and S5B). This shift in gene 
expression profiles suggests that their putative roles may change as 
the larvae shift from the early (1-3 hpi) to later stage (12 hpi) of 
coral morphogenesis. Furthermore, the 1 hpi treatments were 
clearly distinct from those after 2—3 hpi, indicating that genes are 
being regulated earlier than previously reported [6], thus further 
supporting a progression through time with respect to gene 
expression (Figure 2C-D). 

In both assays, principal component analysis and the PCA score 
plots associated with HO org and HO aq (Figure 2E-F) showed a 
clear distinction between the two treatments, even though both 
ultimately resulted in settlement. While the expression profiles of 
12 GOI (mostly immunity-related genes, Figure 1) were found to 
be similar following exposure to the two chemically distinct CCA- 
derived cues, the expression profiles of 15 GOI (Figure 1) were 
different, that is they were regulated following exposure to only 
one of the CCA-derived cues, but were never oppositely regulated. 
The separation of the two treatments was influenced strongly by 
the regulation of T1R-1, AP-1 / cjun, AP-llcFos, MEKK-1, ERK-2, 
Colorless GFP-like, amilGFP and cAMPI ATF6, with a lesser influence 
by JVF-Kb, MAPK-p38, C3, Apextrin, CTL-2 and Tx60-A (Figures 2E- 
F, S4C and S5C). The immunity-related Toll/TLR pathway 
genes contributed mostly to the separation of these treatments. 
These observations raise the possibility that 1) more than one 
component influenced/regulated the immunity-related processes 
involved in both attachment and metamorphosis, or 2) that a 
synergistic effect was causative for this observation. Since the 
cnidarian gene cascade leading to settiement has not yet been fully 
elucidated [10], it cannot be ruled out that numerous pathways 
may be regulated by many different components within the 
environment (i.e. bacterial biofilm, CCA). 

Regulation of the fluorescent protein (FP) genes influenced 
separation between the different time points for JO 10-E exposure 
and also between the two settiement cues. While no initial 
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significant differences were observed 1 hpi with J010-E, there was 
a downward trend in expression levels of amilGFP, amilCFP and 
Colourless GFP-like (Figure 1), with all FP genes down-regulated 
after 12 hpi [12]. FP genes (particularly RFP and CP) were 
reportedly up-regulated in compromised corals, activating an 
antioxidant response [28] and given the strong down-regulation of 
FP genes following exposure to J010-E it is likely the larvae are 
responding to TBP [29] in J010-E. FP genes amilGFP, amilRFP and 
amilCP (Figure 1) were mostly up-regulated 1 hpi with HO aq 
indicating a different role in the early stages of settlement 
compared to metamorphosis. Similarly, Beltran-Ramirez [30] 
reported a strong green FP (GFP) signal in swimming larvae, up- 
regulation of FP in both the ectoderm and endoderm of A millepora 
larvae undergoing settlement and down-regulation in ectodermal 
tissue after metamorphosis [30]. The current study corroborates 
this, providing further evidence that FP genes are differentially 
regulated in swimming larvae post exposure to metamorphosis- or 
settlement-inducing cues and prior to visible metamorphic 
changes. Up-regulation of amilCP 1 hpi with HO aq further 
supports an earlier observation that chromoprotein (CP) gene 
expression was detected primarily during gastrulation and was 
mostly endodermal [30]. Furthermore, A. tenuis larvae exposed to 
CCA chips showed an increase in red FP (RFP) expression prior to 
but not after settlement, while cyan FP (CFP) expression increased 
post settlement [21]. Exposure of A. millepora larvae to J010-E 
resulted in down-regulation of these two FP genes 2-12 hpi which 
correlated to complete metamorphosis. 

The CCA-derived cue HO Qrg represents the organic-soluble, 
hydrophobic, low molecular weight components of the H. onkodes 
extract while HO aq was derived from an autoclaved aqueous 
extract of CCA and contained high molecular weight compounds. 
Both CCA-derived cues ultimately induced complete settlement, 
producing competent polyps that were able to acquire Symbiodinium 
and deposit a skeleton (data not shown). High specificity of large 
molecule-mediated recognition has previously been observed in a 
wide variety of biological systems [3 1] . Based on these observa- 
tions HO aq may be acting similarly to lectins at the cell surface, 
activating the downstream MAPK pathway, whereas HO org 
influences, possibly within the cell, changes in searching behav- 
iour. It should also be noted that like any other surface, CCAs are 
known to harbour bacterial biofilms [32] and that these may 
indeed be the source of the observed inductive capacity of CCA 
fractions on coral larval settlement. Our findings show that two 
different and distinct chemical classes can elicit the same 
settlement response by regulating a different set of genes, 
highlighting that there is more than one pathway involved. 

This study has enabled further elucidation of gene function(s) 
through a time series investigation into their responses to CCA- 
and bacteria-derived cues. Importandy, this gene regulation 
differed not only with regard to the time of exposure but also to 
the morphogenetic outcome; setdement (CCA-derived cues) or 
metamorphosis without attachment (bacteria-derived cue). It 
should be noted that while the genes investigated in this study 
are known for A. millepora, the gene regulation pathways have not 
been fully elucidated and are mostly based on vertebrate systems, 
providing only putative functionality. Even so, the results of this 
study support the notion that regulation of specific genes is driven 
by environmental cues throughout the different phases of A. 
millepora larval settlement, from initial induction to cellular 
differentiation and complete metamorphosis resulting in setde- 
ment. 



Supporting Information 

Figure SI Differential gene expression following expo- 
sure to bacteria- and CCA-derived cues. Significant 
(J)<0.05) change in gene regulation is given as the % difference 
compared to the control. Dark and light green (+) represents those 
genes that were up-regulated (X>20, 20^X<0); orange and red 
(— ) represents those that were down-regulated (0>X& — 20, 
X< — 20). Data normalised to the three most stable genes. Assay 

1 CCA: Rps7, ATF4/5 and MF-kB, Assay 2 CCA: Rps7, Amgalaxin- 
like-1 and Lectin /nemgalB2, Assay 1 J010-E: Rpl9, CTL-2 and 
ATF4/5; Assay 2 J010-E: amilRFP, Rps7 and cAMP/ATF4. Data 
from the 12 hpi experiment (high concentration treatment and 
complete metamorphosis) was taken from Siboni et al. [12], which 
also includes full protein names and description of the genes. 
(TIF) 

Figure S2 Differential gene expression following expo- 
sure to CCA-derived cues, HO org and HO aq , The boxplots 
represent change in gene regulation as % difference from control. 
Only cases which present significant (p<0.05, Kruskal-Wallis 
followed by multiple comparisons of mean ranks) differences from 
control were included. For each gene, the left region represents 
HO org and the right region represents HO aq . 
(TIF) 

Figure S3 Differential gene expression following expo- 
sure to the bacterial cue J010-E over time (1-3 hpi). The 

boxplots represent changes in gene regulation expressed as the % 
difference from the control. Only cases that present significant 
(J)<0.05, Kruskal-Wallis followed by multiple comparisons of 
mean ranks) differences from the control were included. For each 
gene, the left region represents 1 hpi, the middle region represents 

2 hpi and the right region represents 3 hpi. 1 2 hpi data is available 
in Siboni et al. [12]. 

(TIF) 

Figure S4 PCA loading plots of Assay 1 following 
exposure to setdement and metamorphosis cues. PCA 

loading plots (PCI and PC 2) of the % difference in gene 
expression compared to the control for A) J010-E 1-12 hpi and 
HO org/aq , B) J010-E 1-12 hpi and C) HO ()rg/aq . 
(TIF) 

Figure S5 PCA loading plots of Assay 2 following 
exposure to setdement and metamorphosis cues. PCA 

loading plots (PCI and PC 2) of the percentage difference in gene 
expression compared to the control for A) J010-E 1-12 hpi and 
HO org/aq . B) J010-E 1-12 hpi and C) HO org/aq . 
(TIF) 
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